Superconductivity in Ca-intercalated bilayer graphene 
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Recent observation of proximity effect [lj has ignited interest in supercondcutivity in graphene and 
its derivatives. We consider Ca-intercalated graphene bilayer and argue that it is a superconductor, 
and likely with a sizeable T c . We find substantial and suggestive similarities between Ca-intercalated 
bilayer (CeCaCe), and CaC6, an established superconductor with T c = 11.5 K. In particular, the 
nearly free electron band, proven to be instrumental for superconductivity in intercalated graphites, 
does cross the chemical potential in (CeCaCe), despite the twice smaller doping level, satisfying the 
so-called "Cambridge criterion". Calculated properties of zone-center phonons are very similar to 
those of CaC6 . This suggests that the critical temeperature would probably be on the same scale as 
in CaC6- 
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The graphite becomes superconducting after interca- 
lation with alkali elements, with the transition tempera- 
ture ranging from below IK for KCs to 11.5 K for C&Cq 
@; S 0, HI- Mechanism of superconductivity seems to 
be consistent with phonon mediated pairing [6J]. These 
graphite intercalated compounds (GIC) open a promis- 
ing route to an alternative class of superconducting of 
materials with tunable properties. 

The recent discovery of graphene a single sheet of 
carbon atoms, has naturally raised a question of super- 
conductivity in graphene Proximity induced super- 
conductivity in pure graphene has been demonstrated 
[l[. This work shows the potential for developing new 
superconducting devices starting with 2D graphene as a 
basis material. 

The purpose of this Letter is to investigate routes of 
doping graphene with Ca to an extent that would induce 
intrinsic superconductivity. In doing that, we will heavily 
rely on insights from now well understood superconduc- 
tivity in intercalated graphites. 

Obviously, for a large number of graphene layers we 
would have to recover bulk superconductivity of CaCg. 
Hence, if superconducting state is indeed possible, one 
can intercalate multilayer compounds with intermediate 
numbers of layers and investigate the crossover between 
2D superconducting state in a few layers and the increas- 
ingly 3D character of multilayer compositions. This bot- 
tom up approach leads to a reasonable question to ask: 
is an intercalated graphene bilayer a superconductor in 
itself? We will argue that the answer to this question is 
affirmative. 

Let us recall the basics of superconductivity in 
CaC6.@ It was proposed at a very early stage that 
(i) soft Ca modes contribute substantially to electron- 
phonon coupling ll| and (ii) a nearly free-electron three- 
dimensional electronic band, , an analogue of the free 
electron s— band in Ca metal, plays an indispensable role 




Figure 1: Crystal structure of Ca C6 bilayer is shown. 
Ca atom sits in the approximately in the middle of one of 
hexagons of Ca. 



in superconductivityjlOj . An important observation was 
made by Littlewood and collaborators [l(| who pointed 
out that for all known superconducting GIC a nearly 
free electron band (which is well above the Fermi level 
in the pure graphite) crosses the Fermi level (the "Cam- 
bridge criterion" ) . Detailed calculations confirmed both 
conjectures [H, and found that the Caphonons pro- 
vide about half of the coupling strength [l4j , while elec- 
trons experiencing the strongest pairing interaction are 
those in the nearly free electron band, although the car- 
bon electrons also show a sizeable coupling with phonons. 

This naturally suggests Ca as a dopant for graphene. 
We focus on the bilayer graphene case as a first com- 
pound that can be truly intercalated. [Hj]. Moreover, 
intercalation with Ca should provide additional rigidity, 
making the new superconductor structurally robust. It 
is not obvious, however, that intercalating a bilayer will 
be as effective as intercalating the bulk graphite. As 
a very minimum, intercalating a bilayer to the same de- 
gree as graphite provides twice less carriers, e.g., CgCaCg 
vs. CaCg. 16]. Therefore, the first question to ask is 
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whether the emprirical "Cambridge rule" holds, that is, 
whether the Ca-derived nearly-free-electron band does 
cross the Fermi level? Note that as opposed to the 
graphites, where this band is three-dimensional and its 
density of states depends on filling, in bilayer graphcnc 
it is 2D and therefore has filling-independent density of 
states (as long as it crosses the Fermi level). The other 
point is to investigate how similar are the elastic prop- 
erties of Ca intercalated between two graphene sheets to 
that of the Ca in graphite. In what follows we will answer 
both questions affirmatively, based on first principle cal- 
culations, and will therefore suggest that Ca-intercalated 
bilayer graphene should be a superconductor with a crit- 
ical temperature comparable to that of CaCg. 

To this end, we have used the standard full-potential 
linear augmented wave method for band structure 
calculations [TtJ in conjunction with the density func- 
tional theory in a generalized gradient approximation. 
The setup and technical details have been described 
elsewhere 11 1. To imitate an isolated CgCaCg trilayer 
we used an 18.2 A thick slab and optimized the distance 
between the C layers. The latter came out to be 4.64 
A, only slightly expanded compared to that in the bulk 
CaC6, namely. This, by itself, indicates that bonding is 
similar and the phonon property will probably also be 
similar. To verify that latter conjecture we have per- 
formed frozen phonon calculations of a Ca E\ u phonon 
mode, corresponding to the Ca layer sliding with respect 
to the C layers. The results are shown in Fig. [2l 
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Figure 2: Frozen-phonon energy for a E\ u mode calculated 
by displacing Ca along the x direction. The inset shows calcu- 
lated forces acting upon Ca for the same displacements. The 
lines are, respectively, the least square quadratic and linear 
fits. The coefficients agree within 6% and yield for the phonon 
frequences of 120 (123) cm' 1 . (Color online) 

Note that the mode is very harmonic. The calculated 
frequency is 123 cm -1 from the total energy fit, and 120 
cm -1 from the forces, indicating a good convergence. 
This frequency is similar to the corresponding frequency 




Figure 3: Band structure of the bilayer graphene CaeCaCe- 
The size of the symbols indicates the relative C-7r character 
of the electronic states. The arrow "A" points to interlayer 
nearly-free-electron C, band (see text) and the arrow "B" to 
the bonding a band. (Color online). 



in the bulk CaC6|l2|, ss 115 cm , supporting our con- 
jecture that the phonon properties of a Ca layer sand- 
wiched between two C layers are very similar to those 
of Ca in CaCq. This result was computed for purely Ca 
displacements; allowing for the weak hybridization with 
C mode of the same symmetry would bring the frequ ency 
slightly down, making it even closer to that in CaC6 [l2| • 
We have also calculated the two A\ g modes: the in-plane 
C mode, that is expectedly hard, 1370 cm -1 , and the 
one that corresponds to breathing between the C planes. 
The latter mode is absent in the bulk compound. It ap- 
pears to be very soft, 108 cmT 1 . Given that the position 
of the free-electron band strongly depends on the inter- 
layer distance, we assume it will interact strongly with 
that band and thus open an additional electron-phonon 
coupling channel, absent in CaCg. 

Let us now turn to the electronic properties. In Fig. 
[3] we show the calculated bands with C p z character em- 
phasized (the so-called tt bands). Obviously, the Dirac 
points are now well below the Fermi level. More impor- 
tantly, we find the Ca-derived nearly free electron band 
crossing the Fermi level with its bottom located at the 
r point 0.5 eV below the Fermi level. If we place the 
intercalated bilayer graphene on the "Cambridge plot" 
from Ref. :7o[ we observe it to be located above CaC6, 
partially due to the nominally twice smaller doping, 1/6 
vs. 1/3, see Fig. gl 

Other characteristics of the calculated electronic struc- 
ture are shown in Fig. [5] (density of states) and [6] (Fermi 
surface). Several observations are in order. First, the 
density of states per carbon at the Fermi level is very 
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Figure 4: Correlation between the position of nearly free 
band, brought down to chemical potential due to Ca and 
superconducting properties of different materials, from [Icil ]. 
We added to this figure a point (green circle) representing 
C 6 CaC 6 (cf. Fig. ©.(Color Online). 



Figure 5: Density of states projected onto individual atomic 
functions. Note the sharp onset of the bonding C a bands 
(long-dash green line) below —4 eV and the onset of the £ 
band (short-dash blue line) above —0.5 eV. The inset shows 
the total DOS, which is very close to that of the bulk CaC6 
if compared on the per-C basis. (Color online). 



similar to that of the bulk CaC@ : 2.5 states per C12 and 
compared to 1.5 per Cq in CaCfijllj. Second, compared 
to CaCg (see Ref. [ll|, Fig. [6]), the Fermi level appears 
in a minimum of the density of states between two large 
peaks at ~ ±0.5 eV, a feature favorable for crystal sta- 
bility. Finally, it is worth noting that the Ca-projected 
DOS is at least half of the C-projected one. The former 
comes predominantly from the interlayer band. Since 
the LAPW program projects the DOS onto the muffin- 
tin spheres, and the interstitial space in this structure is 
huge, this indicates that the DOS of the interlayer band 
is comparable with that of the C-7T bands, despite the 
fact that the corresponding Fermi surface (small orange 
circle around the T point in Fig. [6J is so small compared 
to the Fermi surfaces of CaC6 and YbCe (Refs. [rl [lij]). 
This is a manifestation of the two-dimensional character 
of the band structure and independence of the DOS (and 
the corresponding electron-phonon coupling [l9j]) of the 
doping level. A useful, albeit elementary, excersise is to 
compare DOS of a 3D parabolic band at a given filling 
with that of a 2D parabolic band with the same mass. It 
is easy to see that the latter is larger as long as the Fermi 
vector in the 3D band kp < tt/c, where c is the interlayer 
spacing in the 3D system (geometrical meaning: as long 
as the Fermi sphere diameter is smaller than the interpla- 
nar distance in the reciprocal space). Note that in CaC6 
and YbC6 this condition is not satisfied, and therefore 
going from 3D to 2D is beneficial for superconductivity. 

As discussed, the coupling of the interlayer band with 
phonons is essentially doping-independent (as shown in 
Ref. [l9j], this holds even when the FS is so small and 
the Kohn anomaly so strong that the phonon self-energy 




Figure 6: 2D Fermi surface of CeCaCe- Note a perfect circle 
around the V point (orange online), derived from the ( band, 
and the triangular Fermi contours derived from the Dirac elec- 
trons. (Color online). 



should be calculated self-consistently including the feed- 
back effects). Based on this, we suggest that pairing 
coupling constant in CagCaCg is similar to the parent 
3D compound. However, this does not mean that ad- 
ditional doping does not help in terms of enhancing the 
coupling. As we know from the CaCg calculations, about 
half of the total coupling comes from the C ir bands, and 
these, being strongly non-parabolic (a triangular shape of 
the largest Fermi surface in Fig. [5] attests to the fact that 
the C bands are still fairly close to the Dirac dispersion) , 
do show energy-dependent DOS, and, by implication, 
energy-dependent electron phonon coupling. Adding sur- 
face Ca atoms to the intercalated bilayer graphene has 
potential to further increase the coupling and enhance 
the critical temperature. 
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Our observation about the nearly free electron band 
crossing the Fermi level (Figs. I3I4| and the similarity 
of the phonon spectra in the intercalated bilayer and in 
the bulk CaCg are the main results of this Letter. Based 
on these similarities, we conjecture that Ca intercalated 
graphene bilayer is a superconductor. 

A few comments are in place regarding the practicality 
of these materials. We know now that the single and the 
bilayer graphenes are electronically inhomogeneous, as 
was seen by scanning probes 2(3, 21 1. Effects of charge 
inomogeneity on the superconducting state would need 
to be addressed in details if indeed superconductivity is 
observed in these materials. Variations of local charge 
density in bilayer graphene would play a role of nonmag- 
netic impurities. As such, they would be subject to the 
Anderson theorem and thus would not be pair-breaking 
for the isotropic or nearly-isotropic s-wave superconduc- 
tivity, as it is believed to be the case in CaC6[L3| and 
presumably is in CgCaCg. Another important effect left 
outside of the scope of this letter is the role of substrate. 
We assumed that substrate effects would be the strongest 
for a single layer. One would need to model the substrate 
in a manner consistent with the DFT approach used here 
to address its role specifically. 

In conclusion, we propose to search for superconduc- 
tivity in a Ca intercalated graphene bilayer. Our esti- 
mates of the phonon frequencies and possible electron- 
phonon coupling constant make this material a plausi- 
ble candidate for superconductivity with T c in the range 
of few K, and possibly above 10 K. We start with an 
undoped bilayer graphene that is nonsupcrconducting. 
Ca intercalation renders a band structure that is rather 
close to that of the bulk CaC6, a known superconductor 
with T c = 11.5X.0, H i, Frequencies of zone cen- 
ter phonons are very close to those in CaC6. The nearly 
free electron band crosses the Fermi level in CeCaCg, 
just as it does in all known superconducting interca- 
lated graphites [3], thus satisfying the "Cambridge con- 
jecture", Fig([4]). Despite the lower doping level, the 2D 
character of this band provides an even higher density of 
states than in CaC6- 

The path to design superconducting material proposed 
here is a bottom up approach, similar to the multilayer 
superconducting films. If the ideas presented here turn 
out to be relevant we might see a new approach to up- 
warard scaling superconducting materials. More broadly 
these ideas fall within the approach of materials by de- 
sign, where we are trying to design materials that target 
particular function, in this case superconducting proper- 
ties with highest T c . 
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